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Abstract 

While  most  previous  work  in  planning  manipulation  tasks  has  relied  on  the 
assumption  of  quasi-static  conditions,  there  can  be  situations  where  the  quasi- 
static  assumption  may  not  hold,  and  the  assumptions  about  the  environment 
must  be  relaxed.  This  is  true,  for  ex^ufnple,  in  a  situation  where  objects  are 
making  and  breaking  contact  at  high  enough  velocities  that  contact  dynamics 
play  a  significant  effect  in  the  motion  of  the  colliding  objects. 

There  has  been  some  work  studying  models  of  collision,  in  particul2u  for  the 
design  and  analysis  of  systems  with  intermittent  constraints,  and  for  the  design 
of  juggling  robots.  Our  work  extends  previous  studies  in  planar  juggling  to  the 
case  of  a  polygonal  object,  using  the  model  of  rigid  body  impulsive  collision. 
Simulations  verify  the  results  of  a  linearized  anedysis. 
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1  Motivation 

Most  previous  strategies  for  plsumiug  manipulation  tasks  have  relied  on  an  as¬ 
sumption  of  quasi-static  mechanics  in  the  analysis  of  the  physical  system.  This 
constrains  the  plans  to  situations  that  are  slow  moving,  and  in  which  contact 
dynamics  can  be  neglected. 

One  can  imagine  situations  these  assumptions  do  not  hold  or  when  a  model 
of  the  contact  dynamics  would  be  useful.  In  catching  an  already  moving  or 
accelerating  object,  for  instance,  the  inertial  properties  of  the  object  affect  the 
motion  which  results  from  the  applied  forces  of  collision.  Knowledge  of  the 
magnitudes  as  well  as  the  direction  of  forces  and  velocities  becomes  important. 
Juggling  and  table  tennis  are  two  such  domains  that  have  been  explored  in 
robotics.  Catching  of  tossed  objects  is  a  related  task  in  which  such  knowledge 
is  useful. 

Another  domain  in  which  such  a  model  may  be  useful  is  in  the  manipulation 
of  objects  by  sliding  on  a  frictional  support  surface.  Much  work  has  been  done 
in  the  analysis  of  quasistatic  pushing  in  the  presence  of  friction  ([13],  [12],  [14]). 
[8]  shows  that  the  motion  of  an  object  on  a  frictional  support  surface  can  be 
determined  if  the  pressure  distribution  of  the  object  is  known.  [12]  and  [14] 
analyze  this  situation  when  the  pressure  distribution  is  not  known.  The  2malysis 
of  [8]  implies  that  for  large  enough  applied  forces,  the  motion  of  the  object  is 
essentially  given  by  the  acceleration  due  to  the  applied  forces.  Since  the  impact 
model  assumes  that  at  the  moment  of  collision  the  impact  force  “swamps  out” 
all  other  forces,  the  use  of  controlled  collision  can  be  useful  in  situations  where 
the  pressure  distribution  is  not  completely  known,  particularly  if  the  support 
friction  is  fairly  low. 

Some  work  has  been  done  to  study  models  of  dyn2Uiiic  collision  for  use  in 
robotic  domains.  [15]  designed  a  dynamically  stable  hopping  robot,  modelling 
the  bounce  as  a  spring  and  damper  system  with  perfectly  inelastic  collision. 
This  system  was  further  analyzed  by  [9].  [1]  designed  a  ping-pong  playing  robot 
which  used  a  simple  model  of  point-mass  collision  to  predict  the  motion  of  the 
badl  after  striking.  [19]  attempts  to  characterize  the  queditative  behavior  change 
in  the  motion  of  objects  upon  collision.  [20]  simulates  and  analyzes  systems 
with  intermittent  constraints,  and  uses  models  of  those  systems  in  planning 
manipulation  tasks.  [5]  analyses  and  designs  a  planar  puck  juggling  system. 
[16]  extends  this  to  the  3-D  case.  This  work  continues  those  studies,  extending 
the  planar  puck  juggling  work  of  [5]  to  objects  with  extent  and  orientation.  The 
intent  is  to  evaluate  the  utility  of  the  impact  model  for  manipulation  planning. 


2  The  Problem 

We  have  a  planar  object  on  a  frictionless  inclined  plane,  pulled  by  the  influence 
of  gravity  down  to  a  movable  “table” ,  against  which  it  bounces  with  coefficient 
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Figure  1:  Basic  System  Geometry 

of  restitution  e.  The  coefficient  of  friction  between  table  and  object  is  /it.  The 
contact  between  the  table  and  the  object  is  a  point  contact,  at  a  known  point 
on  the  object.  We  assume  an  impulsive  impact  model  as  described  by  [17]  and 
used  by  [19],  [20]. 

The  object  is  parameterized  by  (x,y,pd),  where  (x,y)  are  the  coordinates 
of  the  center  of  gravity  of  the  object,  6  is  the  orientation  of  the  object,  and 
p  is  the  radius  of  gyration  of  the  object  ([6]).  The  desired  orientation  of  the 
object  will  be  set  to  =  0.  For  simplicity,  we  can  abstract  away  the  dimensions 
of  the  actual  object,  and  think  of  the  object  as  a  rod,  whose  center  of  gravity 
is  the  center  of  gravity  of  the  object,  located  at  length  /  from  the  from  the 
contact  point^  (See  figure  1).  We  want  to  have  the  object  bouncing  to  a  fixed 
height  while  maintaining  the  desired  orientation.  Ideally,  we  would  prefer  the 
bouncing  to  occur  in  a  fixed  (impact)  position  in  the  ji  j/i  plane,  but  we  will 
initially  disregard  this  constraint,  and  study  the  simpler,  lower  dimensional 
unconstrained  case. 

In  this  examination,  we  will  also  assume  perfect  sensing  and  perfect  control 
of  the  motion  of  the  table,  for  the  purpose  of  examining  the  question  of  whether 
the  desired  behavior  is  achievable  in  theory,  before  exploring  problems  in  actual 
implementation. 

2.1  The  Lossless  Case 

We  first  look  at  the  simplest  case,  the  case  e  =  l,p  =  0.  In  this  case,  the 
equation  describing  the  change  in  velocity  due  to  an  impulsive  collision  can  be 

'Thi*  of  course  ignores  the  question  of  whether  or  not  we  can  actually  strike  the  object 
at  the  angle  desired  at  exactly  the  point  desired,  so  in  effect  we  are  assuming  the  object  is 
“pointy  enough”  at  the  contact  point. 
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written  in  terms  of  the  preimpact  velocity  v~  as 


Av  =  — 2(n^v  )n 


(see  Appendix  A),  where 


n  =  .  COSO  .  (2) 

+  Psin-;3  [  isin/5  J 

As  shown  in  Figure  1,  o  is  the  angle  that  the  table  makes  with  the  horizontal 
in  the  counterclockwise  direction,  /  is  the  length  from  the  center  of  gravity  to 
the  contact  point,  and  0  is  the  angle  that  /  makes  with  the  table  normal  in  the 
counterclockwise  direction.  Note  that  0  =  6  —  a. 

Equation  (1)  simply  says  that  in  configuration  space,  the  object  upon  im¬ 
pact  will  reverse  its  normal  velocity  component  while  its  tangential  component 
remains  unchanged.  Define  the  constraint  surface  to  be  the  set  of  configurations 
for  which  the  object  touches  the  table  without  penetration  for  a  given  table  ori¬ 
entation,  a  (see  [11]  for  a  discussion  of  configuration  speice).  Then  according  to 
equation  (1),  the  object  will  reflect  about  the  normal  to  the  constraint  surface. 
This  is  a  fairly  intuitive  extension  of  the  usual  example  of  a  perfectly  elastic, 
frictionless  point  mass  impact  against  a  flat  barrier. 

In  order  to  get  some  idea  of  possible  solutions  to  our  original  problem,  we 
can  first  look  at  the  case  of  this  point  mass  bouncing  against  our  table.  If  the 
ball  makes  impact  with  a  horizontal  (o  =  0)  table  with  a  velocity  vector  at  angle 
0  =  9  with  the  vertical,  it  will  leave  after  impact  with  angle  —0,  fall  under  the 
influence  of  gravity,  and  then  (if  the  table  remains  at  the  same  height)  strike  the 
table  again  at  angle  9.  We  would  like  to  control  the  strikes  so  that  9  eventually 
goes  to  zero,  and  the  ball  bounces  straight  up  and  down.  One  way  to  do  this  is 
tilt  the  table  while  simultaneously  moving  it  so  that  impact  always  occurs  at  a 
given  height  in  the  vertical  plane,  say  y  =  0. 

Suppose,  as  in  figure  3,  that  the  ball  first  makes  impact  with  velocity  at 
angle  9~  to  the  vertical.  If  the  table  is  tilted  at  angle  a  at  impact,  the  ball  will 
leave  the  table  at  angle  —0  =  —{9~  —  q),  which  is  equivaleti  o  angle  —9~  +  2a 
to  the  vertical.  Hence,  on  the  next  impact,  it  will  strike  with  the  negative  of 
that  angle  with  respect  to  the  vertical.  This  gives  the  system  equation 

^n  +  l  —  ~  2Q:n  (3) 

where  9„  gives  the  velocity  angle  with  respect  to  the  vertical  just  before  the  the 
nth  impact.  If  we  choose  a  to  be  proportional  to  9,  say  o  =  k9,  the  recurrence 
relation  becomes 

<?„+!  =  (1  -  2k)9„  (4) 

or 

^„  =  (1  -  2K)'*f?o.  (5) 
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Figure  2;  Bouncing  Badl  example 


This  will  eventually  drive  the  system  to  its  equilibrium  state,  6  =  0,  as  long 
as  |1  —  2k|  <  1,  or0<«<  1.  We  can  try  to  extend  this  strategy  to  the  rigid 
body  case. 

2.1.1  System  Equations 

We  can  rewrite  equation  1  in  the  form 

v+  =  (I  -  2nn^)v-  (6) 

where  V  is  the  velocity  vector  immediately  before  impact,  is  the  velocity 
immediately  after  impact,  and  I  is  the  identity  matrix.  After  impact,  during 
the  ballistic  phase,  we  have  the  equations 

x(t)  =  Xo 

y(<)  =  yo-Gt  (7) 

pe{t)  =  pBo 

where  the  subscripted  velocities  Me  the  velocities  at  the  beginning  of  the  ballistic 
phase,  ».  e.  the  velocities  given  by  v"*" .  G  is  the  acceleration  of  gravity.  The  next 
impact  occurs  at  j/contoc»  =  0,  when  the  contact  condition 

/  cos  ^0  +  yot  -  ^Gt^  =  /  cos(0o  +  ^oO  (8) 

is  satisfied.  The  time  until  next  impMt  is  a  function  of  ^o,  l/o.  which 

in  turn  are  functions  of  the  configuration  just  prior  to  impact,  (S",v~).  Let 
the  configuration  at  impact  n  be  given  by  x„  =  {p$,  x,  y,  p6)„,  and  call  the  time 
of  next  impact  t(x„).  Then  the  impact  equations  plus  the  contact  condition 
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Figure  3;  BaJl  bouncing  against  tilted  table 
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describe  a  nonlinear  recurrence  relation. 


x„+i  =f(a,x„,r(x„)).  (9) 

If  we  set  a  =  then  f  is  completely  a  function  of  x„,  and  we  can  try 
to  find  a  fixed  point  x*.  By  inspection  of  (6),  (7),  and  (8),  we  can  see  that 
($  =  $  =  0,y  =  y* =  X*)  defines  a  set  of  fixed  points,  with  r  =  —2y‘fG. 
Further,  we  know  that  for  a  lossless  system  that  the  energy  at  contau:t,  + 

y^+p^6^)+mGl  cos  6,  is  constant  from  impact  to  impact,  so  our  initi2d  conditions 
determine  the  energy  surface  to  which  we  are  constrained.  For  specificity,  we 
choose  an  equilibrium  point  from  our  set  of  possible  fixed  points  to  study.  Since 
we  would  like  our  object  to  bounce  in  place  horizontally,  we  choose  i*  =  0. 
Then,  if  we  drop  the  object  from  a  (center  of  gravity)  height  yo,  with  a  initial 
velocity  vector  vq,  the  corresponding  value  for  y‘  is  given  by 

y*  =  -^2G(yo-lcos<l)  +  vo’’vo.  (10) 

We  wish  to  determine  the  stability  of  the  state  x*  =  (0, 0,  y‘,  0);  if  we  start  off, 
not  exactly  at  this  point,  but  merely  near  it,  can  we  balance  the  object? 

2.1.2  Linearized  Analysis 

Given  x*  =  f(x*)  for  the  recurrence  relation  Xn+i  =  f(x„),  Taylor  expansion 
about  X*  gives  us 


{(x'  +  Sx„)  ~  x’-hSx„+i 

=  x’ +  J(x‘)Sx„  +  h.o.t  (11) 

where  J"(x)  is  the  Jacobian  of  f()  evaluated  at  x.  Ignoring  higher  order  terms 
gives  the  approximate  linear  system 

#x„+i  w  J'(x*)(Jx„,  (12) 


6x„  =  [J(x*)r5xo.  (13) 

From  the  linearized  systems,  we  can  try  to  predict  certain  properties  of  the 
nonlinear  system  in  the  neighborhood  of  x*  (see  Appendix  B). 

The  linearized  system  around  x*  is  given  by  the  system  matrix 


J  = 


0  0 

^  10  0 

0  0  10 

0  0  1 


(14) 
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which  has  eigenvzdues  A,  given  by 


1, 

1, 


Gf^ 


± 


STi - 


(15) 


If  all  the  eigenvectors  Vj  are  distinct,  then  the  system  solution  is  given  by 


4 

6x„=5]c.ViA<",  (16) 

t=i 

the  a  being  functions  of  the  initial  conditions.  Clearly,  in  order  for  the  solution 
to  be  stable,  all  the  eigenvalues  must  be  contained  in  the  closed  unit  disc  of  the 
complex  plane.  Analyzing  A3, 4  for  different  values  of  k  gives 


«  <  1  : 

|A3|>1.|A4|<1 

Unstable 

(17) 

Gp^ 

l•^3.4|  =  1 

Stable 

(18) 

«  >  1  +  . 

|A3|<1,|A4|>1 

Unstable. 

(19) 

The  boundaries  of  the  region  given  by  (18)  are  special  cases,  because  for 
those  values  of  /c  not  all  the  eigenvectors  are  distinct,  and  the  amalysis  is  more 
complicated  (see  Appendix  C).  For  the  region  of  stable  k,  (16)  can  be  written 
in  the  form 


<5x„  =  Cl 


+  C2 


+  2Re 


(20) 


and  the  initial  conditions  give  for  the  constants  of  proportionality; 


Cl 

II 

0 

(21) 

C2 

=  Syo, 

(22) 

C3 

p69q  ..ap69o 
~  2h  2bh 

p69o. 

2b  ’’ 

(23) 

The  linear  an2tlysis  predicts  that  in  the  neighborhood  about  the  equilibrium 
point,  if  there  is  any  deviation  from  y*,  it  will  stay  constant;  if  =  0  and 
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00  [rad 

00  [rad/s] 

yo  [m; 

K 

^critical 

Case  1 

0.3 

0 

0.2 

1.1 

1.5 

Case  2 

0 

0.1 

0.2 

1.1 

1.5 

Case  3 

0.1 

0 

0.2 

1.55 

1.5 

Table  1:  Cases  simulated  for  Lossless  System,  zq  =  0  and  y  =  0  when  object 
was  dropped. 

66q  =  0,  then  6x,S6,  and  6$  will  all  oscillate  about  the  origin  at  frequency  a, 
with  amplitudes  determined  by  S0o-  These,  then,  are  the  initial  conditions  that 
determine  the  stability  of  x*.  If  either  6xo  or  SSo  &re  nonzero,  there  will  be  a 
net  X  velocity,  and  the  object  will  remain  balanced,  but  travel  horizontally  as 
it  bounces.  Strictly  speaking,  this  is  not  really  stable,  since  z  (which  we  have 
been  ignoring  up  until  now)  can  increase  without  bound.  But  since  x*  does  not 
contain  z,  it  does  remain  bounded.  For  the  case  k  =  Kcriueai  =  l+(Gp*)/(y*^0> 
linear  analysis  predicts  instability  of  the  system  (see  Appendix  C). 

2.1.3  Empirical  Verification 

Simulation  showed  that  the  linear  approximation  predicted  the  stability  of  the 
system  reasonably  well  for  different  choices  of  k:  the  region  described  by  (18) 
was  indeed  stable  for  small  initial  velocities  and  for  angular  deviations  up  to 
about  dbO.4  radians  23*).  When  the  initial  impact  angle  was  small,  about 
±0.15  radians  (%  8.6*)  or  less,  the  system  was  stable  all  the  way  up  to  and 
including  k  =  KeriUeai-  When  the  initial  impact  angle  was  about  in  the  range 
±(0.15  to  0.4)  radians,  k  had  to  be  much  closer  to  unity  for  stability.  The  k  =  1 
case  is  always  unstable,  but  the  k  =  KeriUeai  case  can  be  stable,  despite  the 
linear  prediction.  Note  for  comparison  that  the  range  ±0.5  radians  (ss  30*),  is 
the  range  over  which  the  linear  approximation  sin  z  ft;  z  holds. 

All  examples  shown  are  for  simulations  of  an  isosceles  triangle  of  uniform 
mass  distribution,  10  centimeters  wide  at  the  base,  and  20  centimeters  high. 
The  radius  of  gyration  about  the  center  of  gravity  for  this  triangle  is  about 
9.428  centimeters. 

The  figures  for  cases  1  and  2,  which  were  stable,  show  projections  of  the 
system  orbit  in  the  <  9  9  >  plane  {6  is  on  the  horizontal  axis)  and  the  projection 
in  the  <  x  y  >  plane.  Notice  in  Figures  5  and  6  that  the  <  0  6  >  projection 
is  centered  about  the  origin  in  both  cases,  showing  that  0  and  0  oscillate  about 
zero.  In  case  1  (Figure  5),  where  there  were  no  initial  velocities  (when  the 
object  was  dropped),  z  is  also  centered  about  the  origin,  and  y  is  bounded  in 
the  neighborhood  of  yo,  the  y- velocity  at  first  impact.  In  case  2  (Figure  6), 
where  there  was  some  initial  angular  rotation,  z  is  no  longer  centered  about 
zero.  Case  3  (Figure  7),  where  k  >  fteriueai,  was  unstable,  and  the  curve  in  the 
<  0  0  >  plane  eventually  goes  unbounded. 


I 
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Figure  6:  Lossless  case  2  with  =  Orad,So  =  0.2rad/s,  k  =  1.1 


2.2  Extensions  to  More  General  Cases 

We  can  extend  the  above  analysis  in  order  to  generalize  the  lossless  case. 


2.2.1  Inelastic  collision 


For  the  case  e  1  we  can  simulate  the  lossless  case  by  giving  velocity  to  the 
table;  if  we  assume  that  masstabie  >  >7iasso»;eet,  then 

v+  =  (I  -  (1  +  c)nn’')v“  +  (1  +  c)nn^Vfa»je  (24) 

(see  Appendix  A)  and  vtabu  remains  unchanged  due  to  our  assumption  about 
the  relative  masses.  If  we  set 

1  —  e  _ 

VtoM«  =  ■  (25) 

1  +  e 

then  (24)  reduces  to  the  equation  for  the  lossless  case.  This  choice  of  vtabu  cw 
be  thought  of  as  having  the  table  tracking  the  orientation  of  the  object  and  the 
velocity  of  the  contact  point,  and  striking  the  object  with  the  table  oriented  at 
angle  a,  and  moving  in  its  normal  direction  at  the  appropriate  speed.  This  czm 
be  seen  by  noting  that 

ieontocl  =  i  + Id  cos  9 

y contact  =  y  +  l9sm9,  (26) 


when  measured  in  the  global  frame  (corresponding  to  a  =  0).  Then  another 
choice  for  Vtat/e  that  will  give  the  same  system  equations  as  (25)  is 


1  “  S  «  .  T 


^contact 

ycontact 

0 


(27) 
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Figure  7;  Unstable  case:  k  >  Kcriucai 
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Bo  [rad]  6q  [rad/s] 

yo  H  y’  M 

y'  [m/s] 

K 

^critical 

ke 

Edesired  [Nni] 

0.3 

0 

0.2 

0.25 

-1.27 

1.1 

1.28 

0.165 

1.73 

0 

0.1 

0.2 

0.25 

-1.27 

1.1 

1.28 

0.165 

1.73 

Table  2:  Cases  simulated  for  system  with  energy  feedback,  e  =  0.5 


2.2.2  Energy  Feedback 

If  in  addition  to  having  the  object  bounce  straight  up  and  down,  one  also  wanted 
the  object  to  bounce  to  a  specific  (center  of  gravity)  height,  y*,  the  associated 
(unit  mass)  energy  level  r;*  =  Gy'  C2ui  be  used  as  an  additioned  feedback  term. 
As  in  [5] ,  we  can  use  the  feedback  law 

Stable  =  -  +  «£(•?*  -  7?)j  \ contact, 

n’  =  Gy'  (28) 

=  GI  +  ^y'^,  at  impact 

T]  =  Glcos8+^y^  at  impact. 

For  stability,  ke  should  be  in  the  range  0  <  Kg  <  =  >iEmat  ([5]).  o 

still  is  set  to  k6,  and  the  new  critical  value  for  k  is  now 

«er<.<eaJ  =  1+  (29) 

The  eigenvectors  and  eigenvalues  of  the  linearized  system  are  essentially  the 
same  (for  0  <  k  <  HeriUeai),  except  the  eigenvalue  corresponding  to  the  eigen¬ 
vector  [0  0  1  0]^  (which  is  the  eigenvector  corresponding  to  6y)  is  now  given  by 
the  value  (1  —  tj;^)  <  1>  reflecting  the  linearized  prediction  that  the  devia¬ 
tion  in  y  goes  to  zero,  i.e.  that  the  system  will  converge  to  the  correct  energy 
surface^.  Experiments  confirm  that  for  low  values  of  ke  (about  O.lKsma*)  the 
energy  does  indeed  converge  to  the  correct  level  (as  seen  in  Figures  8  and  9  for 
cases  4  and  3),  and  then  the  behavior  of  the  system  is  similar  to  the  lossless  case. 
For  values  of  ke  much  higher  than  O  lKEmax,  the  system  is  generally  unstable. 

2.2.3  Friction 

If  we  remove  the  assumption  that  fitaHe  =  0,  the  impact  equations  become 
nonlinear,  reflecting  the  nonlinearity  of  Coulomb  friction.  Although  the  analytic 
approach  becomes  more  difficult,  empirical  studies  for  various  values  of  y  found 
this  case  to  be  unstable  when  using  the  table  tilt  rule  explained  above,  somewhat 

^Strictly  (peaking,  because  this  is  no  longer  a  lossless  system,  the  fact  that  three  of  the 
eigenvalues  are  on  the  unit  circle  means  we  cannot  use  the  linearized  system  to  rigorously 
prove  anything  about  the  nonlinear  system  (see  Appendix  B),  although  [5]  does  present  some 
complicated  stability  arguments  for  the  point  mass  case.  We  can  also  use  intuition,  our 
knowledge  about  the  lossless  case,  and  empirical  evidence  to  help  us. 
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Figure  10:  Tangential  and  Normad  Impulses 


contrary  to  expectations.  Apparently,  the  law  cannot  compensate  for  the  energy 
lost  in  the  tangential  direction,  and  in  fact  often  added  an  impulse  in  a  direction 
that  increased  6,  contributing  to  the  tip-over  of  the  object.  This  cam  be  seen  by 
looking  at  the  equation  for  the  moment  due  to  the  impulse  (see  figure  (10)) 

M  =  Pnl  sin  0  +  Ptl  cos  0 

=  Pnl sin 0  +  /iPnSgn  -  Vtel  cos 0  (30) 

for  the  case  of  sliding  contact.  Here  Pt,Pn  are  the  tangential  and  normad  com¬ 
ponents  of  the  contact,  amd  VtC  is  the  tangentiad  contact  velocity.  For  small 
angles,  sin/3  ss  /?  and  cos «  1  —  /?^  »  ;9,  so  when  0  is  smadler  thap  n,  the 
moment  due  to  frictional  forces  cam  potentiadly  cancel  out  the  desired  moment, 
and  cause  the  object  to  rotate  in  the  wrong  direction.  This  problem  can  be  com¬ 
pensated  for  somewhat  by  increasing  k,  which  in  general  causes  0  to  be  lau-ger. 
Increasing  ke  also  prolongs  the  time  that  the  object  cam  be  kept  upright,  but 
the  tangential  forces  increase  the  horizontal  motion  of  the  object,  and  hence  the 
energy  dissipated  to  friction,  and  eventually  the  object  falls  over.  This  difficulty 
cam  probably  be  circumvented  by  striking  the  object  at  an  amgle  that  minimizes 
tamgential  velocity,  or  by  striking  at  a  different  place.  Neither  of  these  options 
has  been  explored  as  yet,  since  both  apparently  require  that  more  attention  be 
paid  to  the  actual  dimensions  of  the  object. 

2.2.4  Imperfect  Sensing  and  Control 

Although  the  system  has  only  been  simulated,  not  actually  built,  an  attempt 
was  made  to  approximate  imperfections  in  sensing  and  control  by  adding  some 
Gaussian  noise  to  the  calculation  of  impact  time  used  by  the  table.  This  changes 
the  angle  and  velocity  of  the  table  at  impact  time,  as  well  as  the  y-height  at 
which  contact  is  made.  Zero  mean  gaussian  noise  with  a  standard  deviation  of  10 
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Figure  11:  Case  1  with  noise 


ms.,  which  is  about  an  order  of  magnitude  less  them  the  time  between  collisions, 
was  used.  Although  the  motion  of  the  object  with  noise  added  appeeus  much 
less  smooth,  the  system  remains  stable  in  a  similar  range  of  initial  conditions 
as  the  noiseless  case.  Cases  1  and  4  were  resimulated,  this  time  with  noise,  and 
shown  in  Figures  11  and  12. 


3  Conclusions 

The  experiments  show  that,  under  ideal  conditions,  control  of  the  object  is 
theoretically  possible  using  knowledge  of  the  collision  parameters.  The  biggest 
problems  are,  of  course,  dealing  with  friction  and  setting  up  the  necessary  sens¬ 
ing  and  control.  Although  indications  eire  that  stability  can  be  maintained  for 
reasonably  accurate  robot  and  sensors,  handling  friction  will  be  necessary  before 
this  scheme  can  be  considered  feasible. 

In  terms  of  applicability  to  other  domains,  the  results  may  also  be  useful 
in  the  dual  problem  of  planning  the  acquisition  of  a  stationary  object  with  an 
accelerating  hand,  or  in  catching.  For  these  tasks,  we  consider  the  desired  stable 
state  to  be  zero  relative  velocity  (and  distance)  between  the  object  and  the  robot 
hand.  Then  we  would  like  to  plan  the  movements  of  the  hand  so  that  the  object 
does  not  fall  or  bounce  away  from  the  hand,  but  instead  eventually  settles  there. 
Work  on  this  is  underway. 

A  Appendix:  Impact  Dynamics 

Following  [17],  the  impulse  equations  for  a  body  of  mass  m  colliding  with  point 
contact  against  a  barrier  (where  it  is  assumed  that  mjamer  >>  m)  are  given 
by 


m(vt  -  vto) 
Tn{v„  -  v^q) 
mp^(u;  —  Wo) 


P, 

Pn 

Pu  —  Ptl  COS  0  -f  Pnt  sin  0 


(31) 


Here,  v,  is  the  relative  velocity  of  the  body  tangential  to  the  contact  normal, 
and  t-n  is  the  component  normal  to  the  contact.  The  subscript  0  designates 
initial  velocity.  Pi  and  Pn  are  the  components  of  the  impulse  in  the  tangential 
and  normal  directions,  respectively.  The  velocity  of  the  contact  point  is  given 
by 

v,c  =  Vt  +  luCOS0 

Vne  =  v„  +  lwsin/?.  (32) 


As  in  Figure  1,  0  is  the  angle  of  the  line  from  the  center  of  gravity  of  the  object 
to  the  contact  point  with  respect  to  the  contact  normal,  in  the  counterclockwise 
direction. 

Combining  (31)  with  (32)  gives  us 


P  cos^  0  sin  0COS0 

Vtc  =  VfcO  -I-  - — ; - Pt  +  - - Pn 


Vne  —  VneO  "I" 


mp*  mp- 

P  sin  0  cos  0  p^  +  P  sin^  0 


mp^ 


-P,+ 


mp^ 


Pn 


(33) 
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Newton  s  mode]  for  impact  (see  [17],  [19],  or  [4])  divides  the  (normal)  impulse 
into  two  parts,  compression,  when  the  colliding  objects  are  moving  into  each 
other,  and  restitution,  when  the  objects  move  away  from  each  other,  Newton’s 
hypothesis  is  that  the  impulse  of  restitution  and  the  impulse  of  compression 
are  in  the  ratio  e.  We  define  P^o  as  the  impulse  at  full  compression:  that  is, 
when  I'nc  =  0,  This  assumption,  plus  the  assumption  of  no  friction  give  us  the 
expressions 


P,  =  0:  P„  =  (l  +  e)P„o. 


(34) 


Substituting  (34)  into  (33)  gives 


Pn 

Pa 


-(1  +  e) 


mp-‘ 


p-  +  P  sin‘  3 


^ncO 


— (1  +  e)/sin;? 


mp 


4-  p  sin'^  3 


;-2 


(35) 


If  now  we  assume  (again  following  Figure  1)  that  the  collision  reference  frame 
is  oriented  at  angle  a  with  respect  to  the  horizontal  (the  x-y  frame),  and  9  is 
the  angle  the  object  makes  with  the  vertical,  then  3  —  9  -  a.  If  we  rotate  (35) 
and  (32)  into  the  x-y  frame,  we  will  get  (using  the  notation  from  equation  (1)) 

Av  =  -(1  +  e)nn^v~,  (36) 


from  which  (1)  follows  in  the  lossless  (e  =  1)  case,  (24)  follows  as  well,  if  we 
recall  that  (36)  refers  to  relative  velocity,  and  that  the  table  velocity  is  escumed 
unaffected  by  the  collision. 


B  Appendix:  Linear  and  Nonlinear  Systems 

Although  the  equations  for  our  system  are  nonlinear,  we  might  hope  to  gain 
some  understanding  of  its  behavior,  at  least  about  x”,  by  the  linear  approxima¬ 
tion,  Nonlinear  systems  theory  ([3])  tells  us  that  if  the  linearized  approximation 
is  asymptotically  stable/unstable  (in  other  words,  if  J  (x* )  has  no  eigenvalues  on 
the  unit  circle),  then  the  nonlinear  system  is  also  asymptotically  stable/unstable. 
This  would  normally  justify  our  examination  of  the  linearized  system.  Unfortu¬ 
nately,  the  system  under  consideration  here  always  has  at  least  two  eigenvalues 
on  the  unit  circle,  so  we  cannot  use  this  standard  argument.  However,  this  sys¬ 
tem  is  conservative,  and  therefore  J  has  determinant  1,  It  can  be  shown  ([10], 
[2])  that  the  following  are  true  in  this  case: 

1,  If  A  is  an  eigenvalue  of  the  system,  then  A*,  the  complex  conjugate,  is  also 
an  eigenvalue^, 

^This  is  eUways  true  for  a  system  with  real  matrix  entries. 
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2.  If  A  is  an  eigenvalue  of  the  system,  then  1/A  is  also  an  eigenvalue. 

Therefore,  for  a  conservative  system,  stability  is  only  possible  if  the  linetur 
system  values  are  all  on  the  unit  circle.  The  case  where  the  eigenvalues  zure  ±1 
is  even  more  problematic  (see  [2]),  but  fortunately  in  the  case  at  heind,  the  two 
complex  exponential  eigenvalues  seem  to  be  the  critical  values  to  examine. 

C  Appendix:  The  case  k  =  Kcriucai 

For  the  case  «  =  Kcriueai,  A3, 4  both  become  1,  and  the  corresponding  eigen¬ 
vectors  both  become  of  the  form  [1/31,^32,0, 1/33],  1/31,1/32,  t'sa  €  3f.  This  means 
that  J  is  no  longer  diagonalizable,  but  contains  more  genertd  Jordan  blocks;  in 
other  words,  the  “generalized  eigenvectors”  must  all  satisfy  either 

JV.  =  A.Vi  or  JV<  =  A.V.  +  V._i  (37) 

([18]).  If  we  find  the  fourth  “eigenvector”  and  expamd  our  solution,  we  obt2un 

0  0  V31  14i 

6x„  =  Cl  Q  +C2  J  -f  (c3(— l)”-l-nc4)  -(-C4(-l)'*  (38) 

0  _  0  _  _  143 .  .  143 

which  would  appear  to  only  be  stable  for  a  limited  number  of  initial  conditions, 
at  best.  However,  this  case  seems  to  be  fzurly  stable  in  simulation. 
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